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Effective Chiral Lagrangian is derived from QCD in the framework of Field Correlator Method. It contains 
the effects of both confinement and chiral symmetry breaking due to a special structure of the resulting quark 
mass operator. It is shown that this Lagrangian describes light pseudoscalar mesons, and Gell-Mann-Oakes- 
Renner relations for pions, eta and K mesons are reproduced. Spectrum of radial excitations of pions and K 
mesons is found and compared to experimentally known masses. 



PACS: 11.15.Tk,11.30.Rd,12.38.-t,12.39.Fe,14.40.Aq 



1. QCD is known to possess two highly nontrivial fea- 
tures at low temperatures, namely confinement and chi- 
ral symmetry breaking (CSB). At some critical tempera- 
ture phase transitions of deconfincment and chiral sym- 
metry restoration occur. From lattice calculations it is 
known that these two phase transitions take place at 
the same temperature 12| ■ The fact that two critical 
temperatures coincide was not fully understood so far. 
*^ This work is a continuation of a series of papers El Eli 
Oh where it is argued that CSB occurs due to confinement 
O 'in a very nontrivial way. 

It was shown in 3 that effective four-quark interac- 
'tion leading to spontaneous chiral symmetry breaking, 
^ .occurs in QCD due to confinement, and is associated 
k> with the QCD string. Thus, CSB is closely connected to 
j_j confinement. In this approach the Effective Chiral La- 
grangian (ECL) containing fields of light pseudo-scalar 
mesons is derived from QCD Lagrangian. This is done 
by integrating out gluon fields and performing boson- 
isation. At the same time confinement is taken into 
account through specific form of gluon-field correlators. 

As a result, expanding in powers of (derivatives) of 
bosonic fields, one obtains the ECL similar to the cele- 
brated Gasser-Leutwyler Lagrangian .6. , however in the 
nonlocal form 

We expand ECL in powers of meson fields, and 
reproduce standard Gell-Mann-Oakes-Renner relations, 
while meson masses are zero in the chiral limit. It 
is shown that the vanishing of meson masses happens 
due to cancellation of two terms in Green's functions 
of mesons. Poles of Green's function corresponding to 
radial excitations of pseudoscalar mesons are displaced 
from the masses, obtained in Hamiltonian approach 
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without CSB effects (see e.g. [7] and references therein), 
and are shifted down by less than 15 %. 

2. We consider Euclidean partition function for quarks 
and gluons in the presence of external classical currents 
Vn, an, s and p 



Z = J DADipDip exp [-(So + Si + Stat + S g . f . + S gh ) 

5 = \j d^x {F^)\ 

51 = —i J d 4 x r ijjf + v + 75a + s + ijsp)f 9 ip 9 , 
Stat = - / d 4 x^ f gA a t a ^ f . 



(1) 

Here /, g — 1,2,3 are flavor indices, t a are generators 
of color SU(3) group, trt a t b = S ab /2, a = 1, ... 8. S g . f . 
and Sgh are gauge fixing and ghost terms. 

Next, we use the generalized contour gauge [SI El 
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Here z v (s, x) belongs to a set of contours, with proper- 
ties: z v (0,x) — xq, z v (\,x) = x u , xq is a fixed point. 
In what follows the exact position of contours is unim- 
portant for our analytical results, while for numerical 
estimates we will assume that contours are chosen to 
minimize the spectrum (and area of the string world 
sheet), to be called the minimal set of contours. 

The reason we use contour gauge is that it allows 
to express gauge field A^ through field strength tensor 
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F^u . Now we are in position to integrate out gluon field 
A^, expressing the result in terms of field correlators: 



(3) 



Z = J J D^exp [-(Si + S eS )} 
exp[-S cS ] = (exp [-Si n t])A 



We use cluster expansion to evaluate this average over 
gluon fields 

(exp [S^U = exp M^i j , 
((Stat)) - (Si*) a = 0, (4) 

((Sint)) = (Sint)A - (5'int)i = (Sf nt )A, 



It is clear that gauge invariant quantities like spectrum 
and Green's functions computed with the help of S e g 
do not depend on the chosen contours, if all terms of 
cluster expansion are retained in Q. In what follows 
we will use the Gaussian approximation, and consider 
only first two terms in cluster expansion, n = 1,2. As 
was shown in |10l II 1| the Gaussian approximation on 
minimal surfaces is accurate within few percent. Thus, 
we have 
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, d 4 xd 4 yip{ a (x)ip f jf3 (x)ip 9 k7 (y)ipf s (y)x 

f 1 d dt dZp ^ S ' X ^ d^(s,3Q dz p/ (t,y) dz y (t,y) 
Jo <9s dx^ dt dy v 

x ([F pA (z( S ,x))] y . [F p , x ,(z(t,y))] kl ) A (r) at3 m lS ■ 

(5) 

Here i,j,k,l are color indices, a, (3,^,6 are spinor indices. 
Inserting parallel transporters xo) and <&(y, xq), 
which are identically equal to unity in contour gauge, 
one finally gets expression for effective action 



S e n = ~ / d A xd 4 y^{ a (x)^(x)^(y)^ s (y)x 



x (SjkSu - -^rSijSki) J a/37S (x,y), 
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N? — 1 



/ dcdt dZp ^ S,X " > dzx ( s ' x } dy( f ' V) dz v {t,y) ^ 



ds 



dx„ 



dt dy v 



x tr (F pX (z(s, x), xq) Fp> X i (z(t, y), x )) A , 
F(u,x ) = §(xo,u)F(u)$(u,x )- 



Performing bosonization, and keeping only scalar- 
isoscalar and pseudoscalar-isovector (corresponding to 



pions, K and rj mesons) terms, one arrives at the quark- 
meson Lagrangian (see 0] for details): 

Z = J Dij)Dil)DM s D(j) a exp [Squ] , 

Sqm =- d i xd i y (x) x 



x (i(d + v + -y 5 a + s + i^pYjpS^ {x-y)+ 



(7) 



+ iM s (x,y)U^(x,y))^(y)- 

-2N f (J(x,y))- 1 M*(x,y) 
J(x,y) = J^(x,y), 

It is now straightforward to integrate out quark fields 
to obtain the Effective Chiral Lagrangian: 



Z = J DM s D4> a e X p [-S ECh ], 



J EGL 



2N f I d 4 xd 4 y(J(x,y))- L M*(x,y) + W ((/>), 



W(<j>) = N c tr\n i{d + v + j 5 a + s + i-f 5 p)+ 

(8) 

Here tr refers to flavor and spinor indices and to space 
coordinates. M s is the effective quark mass opera- 
tor, and 4> a are fields of pseudoscalar mesons (up to 
the dimensional factor 2/f, f is the decay constant, 
4>a = 2n a /f). 

Classical equations of motion are given by 



N B Tr 



N c Ti (-5 (x, 2 /)M s ( a; ,y)e^ t "^^^7 5 i Q ) = 0, 
(iS (t> (x,y)e^ t ^ x ^) - 

W f (J(x,y))- 1 M s (x,y) = 0, 
1 



(9) 



Sj,(x,y) = (x]^ 

This leads to solutions 

^ 0) (z,2/) = 0, 



■|y> 



M^(x,y) 



AN f 



J(x,y)Tr(S(x,y)). 



(10) 



S(x,y) = S^(x,y) 



Second equation in (|1(J|I is a nonlinear equation for Mg , 
and the existence of a nontrivial solution is a mani- 



(6) festation of the chiral symmetry breaking, since M. 



(0) 



is scalar. The system of equations (|10f) for and 
S(x,y) was considered in f° r the special case of 



Pseudoscalar mesons and their radial excitations 



3 



heavy-light mesons, and it was shown that it has a con- 
fining scalar solution for M$ (x, y) ~ cr\x~ x \S^ (x—y) 
for large distance \x — xq\ from quark to (heavy) anti- 
quark at the point xq. It is clear that the same type of 
solution occurs at large interquark distances for light- 
light mesons, which means that confinement and CSB 
occur spontaneously and simultaneously from the non- 
trivial solution of the system (TT77f> . 

3. We consider ECL (JSJ), expanding it in powers of the 
field 4> a up to the second order, and we introduce cur- 
rent quark masses Mf = diag(m„, m d , m s ). Neglecting 
external currents, one obtains 



W(4>) = N c tr\n i(d + M f + Af s e i75t ^) 



=N C tr In i(d + Mf + M s )+ 

-M sl5 t a 4> a - ^M s t a t h (j) a (f)b \ , (11) 

W^(4>) = -^tr[iS-(M s t a t b( j )a( j )b ) + 
+ S ■ {M s t a <j) a ) • 75575 • (M s t b (p b ) 



Taking trace in flavor indices, (see 0] for details) , one ar- 
rives at the following expression for the term quadratic 
in meson fields: 



T^ (2) (0) 



d xd y 



W nw (x,y)(p*(x)(f)- !r (y)- 



+ W K K{x,y)4>* K (x)(t>K{y)+ 
+ W K o K o (x, y)(j)* K o {x)(j) K o (y)- 



(12) 



1 



9 w ^( x ^y)Mx)(t>j(y) 



i,i=3,8 

where, for example, 

W^(x,y) = ^Tr S u (x,y)M s (y)j 5 S d (y,x)j 5 M s (x)+ 
+ iS u (x, x)M s (x)S 4 (x — y) + (u <-> d) . 

(13) 

S u , Sd and S s are quark propagators (|10|l with cur- 
rent mass of corresponding quark in the denominator. 
Here we have taken the local limit of nonlocal opera- 
tors M s (x,y) — > M s (x)S 4 (x — y), 4>(x,y) — > (f>(x), which 
is obtained when gluonic correlation length T g in the 
correlator (FF) tends to zero. 

Two terms in (I13f) correspond to connected and dis- 
connected diagrams, which cancel each other in the zero 
momentum limit. This cancellation is exact in the chi- 



ral limit. To be more precise, quadratic term for zero 
momentum (i.e. when 4>{x) — const) takes the form: 



M/ (2) (0) 
4 



zero momentum 
, : "li i ' ni>d 

d x 



2 

+ 0(m 2 ). 



Tr (-iS u (x, x) - iS d {x, x)) </>* (f> n 



Taking into account that 



1 SZ[v, a, s,p] 
Z 5s(x) 



(14) 



N c Tr(iS(x,x)), 
(15) 



where (ipip)M and (^V')b denote quark condensate in 
Minkovski and Euclidean space respectively, and that 
4> a = 27T a / f, f is the decay constant, ir a are physical 
meson fields, one finds 

fMl ± = 2m\{qq)\ + 0(m 2 ) 
f 2 M 2 = 2m\{qq)\ -e+ 0(e 2 ) + 0(m 2 ) 
f 2 M 2 K± = (m u + m s ) \{qq)\ + 0{m 2 ) 
f 2 M 2 K „ = (m d + m.) \(qq)\ + 0{m 2 ) 
2, 



; {m + 2m s ) \(qq) \ + e + 0(e 2 ) + 0(m 2 ). 



(16) 

Here m = (m u + md)/2. We have neglected differences 
between quark condensates for different flavors, correc- 
tions are of order of m 2 . Small mixing of 03 and (As 
states due to isospin symmetry breaking (proportional 
to m u — md) yields a correction e to pion and rj meson 
masses: 



cos(<5)(/> 3 + sin(<5)</> 8 , 

77s ~ - sm(6)<j> 3 + cos(S)(j) 8 , 

, ^ r- nid — m v 
tan(2<5) = VS- c 



2m s - (m u + md) 
\(qq)\ (m u - m d ) 2 



(17) 



4m s - 2(m u + m d ) 
(5-0.6°. 

Thus, ECL JHJ leads to correct Gell-Mann-Oakes- 
Rcnncr relations for all light pseudoscalar mesons. 

4. Let us now consider Green's functions of mesons, 
generated by the pseudoscalar currents: 



<8 ab (x,y) = (J 5 a (x)Jt(y)) = 

J a( x ) = ^(x)l5ta^(x), 
pf<>(x)= Pa (x)tl°. 



1 



S 2 Z 



Z 5p a (x)5p b (y) : 



(18) 
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From the ECL one obtains 

®ab(x, y) = \\ DM sD(j)a CXp [-SECT,] X 

x N c TT{S^{x,y) l5 t a S^{ y) x) l5 t b )- (19) 

- N 2 Tr (S^(x, x)j 5 t a ) Tr (S^(y, y)^h) 

Taking M s at the stationary point (|10() and expanding 
S$ in terms of <f) around 4>a — 0, one finds: 

S$(x,y) = S(x,y)+ J d 4 zS(x, z)-M s (z)j 5 (j)a(z)t a -S(z,y) 

(20) 

As argued in 0], the coupling constant g-^qq is of the or- 

— 1/2 

der of N c , and thus in large N c limit pion exchanges 
are suppressed. This allows to neglect pion fields in con- 
nected terms (first term in equation l|19|) ). and consider 
only one pion exchange in the disconnected term. Re- 
sulting expression will contain two terms, both of order 
of N c . 

Taking into account, that S(x,y) = 
di&g(S u (x,y),Sd(x,y),S s (x,y)) is diagonal in fla- 
vor, one finds: 

<& n +„+(x,y) = -y Tr (S d (x,y)j 5 S u (y, x)j 5 ) - 
N 2 f , 

~t I rf4z i rf4z 2Tr(S u (x, z 1 )M s {z 1 )^S d {z ll x)^) x 

x Tr (S d (y, Z2)M s {z 2 )^Su{z2,y)lb) G^(z 1 ,z 2 ). 

(21) 

Other Green functions differ only in flavor indices. Here 
G$„(zi,Z2) = {4>rt{ z i)<t>w{ z 2)) is the propagator of pion 
field. This formula can be illustrated with the Feynman 
diagram: 



taken into account that Green function is translational 
invariant, i.e. depends only on (x — y), one has: 



(0,M) 



(24) 



G { °l^{x,y) = Tr (S d (x,y)j 5 S u {y,x)j 5 ) , 

G^+nl ( x ' y) = Tr ( s u(x, v) M s{y)i5S d {y, 37)75) ■ 



d 



®7r+7r+ (X, y) = 



(22) 



To find the pion propagator one should consider 
equations I|12I13() : 

(G* v (x,y))~ 1 = W irv {x,y) = 

= ^Tr[2S d (x,y)M s (y) l5 S u {y,x)M s (x) l5 + (23) 
+ i(S d (x, x)M s (x) + S u {y,y)M s {y))5\x - y) 



Going over to the momentum space <S(x,y) — 
J d i k/(2Tr) 4 exp(ik(x — y))*3(k), where we have also 



Due to Eqs. (f 141 — II fij) the pion propagator (|23|1 has pole 
at k 2 = — M 2 ± and can be rewritten as 

a't> (i,) _ 2l \ 

" +7T+[ ' N c G (0MM) (k) G (0 MM) = _ , 

7T-r7r"r V ' 7T~ t ~7T~ t ~ v TT^ ' 

G^ + M) (x, y) = Tr (S u (x, y)M a (y) l5 S d (y, x)M s {x) l5 ) . 

(25) 

As argued in 4 , all three Green functions G| r + 7r+ , 
have the same set of poles, which 



G 



(0,M) 



and G 



(o,MM) 



are poles of the quark model (i.e. confined qq sys- 
tem without chiral symmetry breaking) in pseudo-scalar 
channel, and can be represented as 



<«U*) = -£ 



n=0 



mi 



„(M) 



n=0 



(26) 



,(M)V 



n=0 



where 



c (M) = /^ M(0Vn(0); 



(27) 



c/? n (r) is the 3D spin-singlet wave function of qq sys- 
tem, and M(0) is a constant related to mass opera- 
tor M s , evaluated in j^j through a = 0.18 GeV 2 and 
T g = 1 GeV" 1 to be M(0) = 148 MeV. Thus one has 
for the pion Green function: 



«*+*+(*) 



tt(Jfc) 



2 (fc 2 + M 2 ± )$(fc)- 



2 / JM) 



*^ E (fc 2 + m 2 )(m 2 -M 2 )' ( 28 ) 
(cl M) ) 2 



$(fc) = £ 



n=0 



{k 2 +m 2 )(m 2 - A/ 2 
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Clearly, the Green function l|28[l has pole at k 2 — 
—M 2 ±, and all poles of quark model are cancelled, since 
the same set of poles appears in functions and 
$(fc). The radial excitations of 7r ± meson are given by 
zeros of the function $(fc). In the first approximation it 
reads: 



k 2 = 
Si = 



1 c 2 (m 2 -m 2 )(m 2 -M 2 ± ) 
m 2 c 2 {ml - M 2 ± ) + c 2 (m 2 - M 2 ± ) 



(29) 



Next, we plug in numbers: 

m u = 0.005 GeV, m d = 0.009 GeV, 
a = 0.18 GeV 2 , 
a s = 0.3, 



= 0.17 GeV, 



(32) 



and taking into account that lowest state is shifted ex- 
actly to physical value of meson mass (due to Gell- 
Mann-Oakes-Renner relations), we finally get the fol- 
lowing chiral shift of reference (quark model) spectra: 



Masses of K°, K° radial excitations can be found 
from H29f) with the exchange of ir meson mass and ref- 
erence spectrum with those for K mesons. Numerical 
results for masses of radial excited states are presented 
in the next section. 

It should be mentioned, that rj meson requires sepa- 
rate consideration, because of it's mixing with isoscalar 
state T)', which is different for mesons and their radial 
excitations. Study of r/ meson spectrum and mixings is 
planned for the future work. 

5. Masses and wave functions of reference spectrum can 
be obtained from the QCD string Hamiltonian (first de- 
rived in ^3 ^] , and improved to take into account 
quark self-energy in JE|), where we have put L = 0: 



H = 



'"-2 
2^2 



Ml + M2 



2fi 



00) 

3 r 



Here mi and m-i are current masses of quarks, \i\ and 
fi2 are einbein parameters, to be found from the eigen- 
values of Hamiltonian l|3UI) via dM n {^\, ^2) / d^i = 0, 
<9M„(/xi,/i 2 )/<9^2 = 0; fx = nin 2 /(pi + M2), and p r is 
the radial component of momentum. This Hamiltonian 
allows to find spin averaged masses and wave functions. 
Spin-spin interaction can than be taken into account as 
a perturbation: 



M n = M n ( Ml)/ i 2 ) + 32 7 sSls V (0)| 2 - 



9/^1^2 



2a ( 1 



M1M2 
1 



(31) 



A SE = 1 )m i)~0.9-rl 

Ml M2 



Ase is the quark self-energy term due to field correla- 
tors. Factor rj is a calculable function of current quark 
masses, but is close to 1 when quark masses are small. 



pions: 

tt(15) 0.51 GeV - 

n(2S) 1.51 GeV - 

tt(35) 2.18 GeV - 

K mesons: 

K(1S) 0.63 GeV - 

K(2S) 1.57 GeV - 

K{3S) 2.21 GeV 



0.14 GeV (exact) 

1.25 GeV (exp : 1.3 GeV) 

1.98 GeV (exp : 1.8 GeV) 



0.49 GeV (exact) 
1.43 GeV (exp : 1.46 GeV) 
-> 2.1 GeV (exp : 1.83 GeV) 



It can be seen, that masses of radial excitations 
are shifted by less than 15%, and the shifts are small 
for high excitations. Moreover, one can estimate that 
SM/M(4S) ~ 0.05 for pions and SM/M(4S) ~ 0.04 for 
K mesons. Also, one can see that masses of higher ex- 
citations and the slope of radial Regge trajectory differ 
from the experimental. The reason is that Hamilto- 
nian H30[) does not take into account effects of string 
breaking, which are important for highly excited states, 
since they have large spatial extent. As was shown 
in J7| the inclusion of string breaking effects does not 
violate the linearity of radial trajectories, which is in 
agreement with experimental data. 3 ) 

6. Effective Chiral Lagrangian JSJ is derived directly 
from QCD Lagrangian in the framework of Vacuum Cor- 
relators Method. This Lagrangian correctly describes 
light pseudoscalar mesons, which are massless in the chi- 
ral limit and satisfy Gell-Mann-Oakes-Renner relations 
when quark masses are nonzero. 

Poles of quark model Green's function are shifted by 
less then 15%, and the shift is small for highly excited 
states. 

The authors are grateful to A.M. Badalian for valu- 
able comments and discussions. 



3 'The notion of linear radial Regge trajectories is connected to 
the Veneziano-type "linear dual models" lljjl . 
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